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Abstract 

We apply the general theory of Cauchy biorthogonal polynomials developed in [2] , [3] to the case associated 
with Laguerre measures. In particular, we obtain explicit formulae in terms of Meijer-G functions for all 
key objects relevant to the study of the corresponding biorthogonal polynomials and the Cauchy two- 
matrix model associated with them. The central theorem proved for such Cauchy two-matrix model is 
that a scaling limit of the correlation functions for eigenvalues near the origin exists and is given by 
certain new random point field, the Meijer-G random field. We conjecture that this random point field 
leads to a novel universality class of random fields parametrized by exponents of Laguerre weights. 
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1 Introduction 

The Cauchy two-matrix model was introduced in [2] as a new kind of an integrable multi matrix model. 
The model consists of a probability measure on the space of pairs Mi , M2 of n x n positive semidefinite 
Hermitean matrices. This probability measure depends on the choice of two scalar functions Vi,V2 '■ 
M+ — > M, called the potentials and is defined as 

1 „-WTr(Vi(Afi) + V-2(Af2)) 

dM(M.,A/.) ^ - dM.dM.. (1-1) 

The parameter iV is a scaling parameter; when considering the asymptotic regime n — > 00 the parameter 
N is also sent to infinity in such a way that — > T S M+. We will assume henceforth T = 1 and that 
N = n. 

It was shown in [2] that, similar to other integrable matrix models, all correlation functions of the 
eigenvalues of the two matrices Mi , M2 can be computed exactly in terms of certain Cauchy biorthogonal 
polynomials (BOPs). The latter consist of two sequences of polynomials {Pn{x),qn{y)}n(^n of exact degree 
n with the defining properties 

Q-N(Vi{x)+V2{y)) 



y 



-pe{x)q,n{y)dx dy ^ Si,n , p„(a;) = c„x" -I- . . . , g„(j/) = c„y" -I- . . . . c„ > 0. (1-2) 



In [3] the algebraic properties of these polynomials were investigated but no concrete example which 
could be considered "classical" was provided. On the other hand, even before the Cauchy BOPs were 
introduced, an instance of such polynomials associated with a classical weight appeared implicitly in [4] 
in the study of the biorthogonal Laguerre ensemble, one of several examples of biorthogonal ensembles 
considered there that allow an explicit computation of correlation functions. In this paper, we apply the 
machinery developed for general Cauchy BOPs in [2], [3] to this particular instance. More precisely, we 
shall consider the case Vi{x) = x + \nx, V2{y) = y + Iny; in this case the scaling factor N in eq. 
(1-2) can be absorbed in a common rescaling of the variables x, y without affecting the general properties 
by setting Vi{x) — (a; -|- a In a;) , V2{y) ~ jj{y + b\ny). The corresponding biorthogonal polynomials 
are thus defined by 

; Pt.{x)qm{y) = ^tm , Pn{x) = CnX"^ + . . . , qn{y) = c^y'^ + . . . . C„ > 0. (1-3) 

x + y 

The present paper has three main goals: 

1. obtain explicit formulae for Pn,qn in eq. (1-3) and for related objects; 

2. find explicit formulae for the correlation functions at finite n; 
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3. formulate a scaling limit of the correlation functions near the origin and thus define a limiting 
random point field. 



We shall show that all the objects naturally appearing in the investigation of the BOPs defined by eq. 
(1-3) are intimately related to Meijer-G functions. It is for this reason that we shall name the limiting 
process mentioned in the last item above the Meijer-G process. 

The results of the paper are summarized in the next section, while the remainder of the paper is 
devoted to the proof of the theorems stated in Section 1.1. 



1.1 Summary of results 

We recall the results of [2] (collected and explained in Appendix B, in particular eq. (B.2b)). The 
correlation functions of the Cauchy two-matrix model are expressed as determinants 



{xi,...,Xk;yi,-.-,ye.) = det 



l<i,j<k 


Hqq {xi^ yj) 

l<i<ka<j<l 


H 11 [yi.xj) 

l<i<l:l<j<k 


H[l\yr,yo) 

l<i,j<i 



(1-4) 



where the kernels -ffp"^ are given by 



4'i\x,y) := x'^y''e-^-yK^\x,v), H^l\x,v) := x'^e-^ K^^\x,v), (l-5a) 

H^^{x,y):^y^e-yK^;^{x,y), H^^\x,y) -.^ K^{l\x,y), (l-5b) 

with the kernels k]jI) defined in terms of the Cauchy biorthogonal polynomials p^, as: 

^ r — y' A ^ 

4:\x,y):=J2p,ix)q,iy), K^^\x,y) -.^ Koi^^v' f \ ,^ , (l-6a) 

y + y 

'Mf^ „^ f T^Mf^' ^^^!e:!^ r.{n). . f ^(n),, ^,,x"^e--' Ax' y">c'y' dy' 1 



K'->{x^y):^ K^->(x',y) W ^ X}r(x,y):= K)->{x',y') 



+ 



x + x' 11 \ 7r2 "'J ' ' x-hx' y + y' x + y 



(l-6b) 



Our first main result is a compact expression for these four kernels for finite n; to present it we set 
a = a + b and define four functions: 

J^2TTir{c+l-~u)r{n + u)r(a-u + l) ' ^ ' 

G^^iw) := /d. r(.-c)r(u)r(n + a+l-.) 
"'''^ ' J 2TTi T{a + l-u)T{n + u) ' ^ ' 
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and 

J^+c 27r« r(l - M)r(c + n + u)T[a - c - u + 1) 
?T / \ f duT{u)T{u + c)Tin + a-c+l-u) ^„ / x /-, mx 

j^+c^m T{a- c+l-u)T{c + n + u) 

In the above expressions 7 is a contour originating at —00 in the lower half-plane and returning to —00 in 
the upper half-plane in such a way as to leave all the poles of the T functions in the numerator containing 
the variable +u inside the contour, while leaving those of the F functions of variable —u outside. Such 
types of integrals are Mellin-Barnes integrals and the expressions above are special cases of Meijer-G 
functions (see [5], 5.3, p. 206). Then 

Theorem 1.1. The kernels H^\k^^ are given by 



Tjin] 
-"00 


\x,y) = 


/ Ha.nitx] 


)HbAty)dt, 


-"01 


\x,y) 


Jo 


-"10 


\x,y) = 


/ Ha.n{tX. 
JQ 


)Hi,,nity)dt, 


^11 


\x,y) 


Jo 


-"■00 


\x,y) = 


Jo 




^01 


\x,y) 


= y'ey 1 




\x,y) = 


x^e^ [ Ga.n(tx)Gb,,,{ty)edt, 
Jo 


Kin] 
"■11 


\x,y) 


= x'^y'^e' 



x + y 



x + y 



The above theorem is the summary of Theorems 3.5,3.6,3.7 that contain the computation of the kernels 
K'fu) , whereas the ensuing expressions for the kernels H^]^ is a simple rewrite using the definitions ( eq. 
(1-5), 1-6, 1-9, and eq. (1-10)). 

1.1.1 The scaling random field at the origin: the Meijer-G random point field 

With the substitutions x — C'^-^^i y — ^.n^^ we arrive in the limit n — > 00 at a novel universality class of 
random point fields that we name the Meijer-G random point field. First we define 

Definition 1.1. Let 

HAC) ■■= CG^iC) - / — ^ ,r", (1-11) 

m := CG.iC) ^ f |^ ^,^^")^^; + '^\ r". (1-12) 
ZiTT i (a — c + 1 — u) 

The contour j is a contour of the form in Fig. 1 linking to all the poles in the numerators of the 
integrands. Note that Hc,Hc are Meijer-G functions as in [5j and 16.17 in DLMF ^ (see definition in 
App. A). 



-'^DLMF=" Digital Library of Mathematical Functions", http://dlmf.nist.gov 
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Theorem 1.2 (Meijer-G random point field and universal class). In the scaling limit the correlations of 
the eigenvalues of Mi, M2 are determined by the two-level random point field with kernels below 

lim n'^H^^ ^\ = goo(C,e), 1™ ^-^iJ^;') \\ = ^oilCO, (1-13) 

lim n-2ij(^) \\ = GioiCO, lim n-'Hi'l^ ^] = GiiiCO, (1-14) 



where 



n—^oo \n n I n— >oo \n n 



Goo{C,0= t HaitCmitOdt, goi{C,0= t HaitOHhitOdt, (1-15) 
Jo Jo 

GioiCO^ f HM)HM)dt GiiiCO^ f HaitQHtitOdt--^ (1-16) 

with He, He as in Definition 1.1. The convergence is uniform for ^,ri within compact sets and the error 
approximation is within 0{n^^). 

Proposition 3.2 provides alternative expressions for the kernels in terms of "point-split bilinear 
concomitants". Section 3.3 is devoted to the proof of Theorem 1.2. We expect the following conjectm'e 
to be true. 

Conjecture 1.1. The Meijer random process obtained in the scaling limit in this paper is universal 
within the class of Cauchy matrix models of the form 

1 p-nTr(yi(Mi) + y2(Af2)) 

A^l{Ml,M2) = — det(Mi)°det(M2)'' — dMidMa 

' ' Z„ ^ ^ ^ ' det(Mi-KM2)" 

with Vi analytic near the origin and the scaling x 1— xn~^. 
1.1.2 Outlook 

The issue as to the real nature of the universality of the Meijer-G random point field for arbitrary 
potentials will be addressed in a forthcoming publication; the proof will be based on the nonlinear 
steepest descent method. An interesting problem which will also be addressed in a separate publication 
is to connect the gap probability of the Meijer-G field to a suitable Painleve transcendent, namely, to 
obtain the analog of the Tracy- Widom distribution of the largest eigenvalue [12]. 

Acknowledgements This paper was initiated at the Banff International Research Station. We thank 
BIRS for the hospitality and for providing excellent work conditions. M.B. wishes to thank SISSA 
(Trieste) for hospitality during which part of the work was completed. M. B. and J.S. acknowledge a 
support by Natural Sciences and Engineering Research Council of Canada. M.G. is partially supported 
by the NSF grant DMS-1101462. 
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2 Setup 



As was observed in [4] the Cauchy biorthogonal polynomials in eq. (1-3) are related to the classical 
Jacobi orthogonal polynomials for the weight x°''^''dx on [0, 1]. We thank A. Borodin for pointing out 
this connection whose main point is as follows: consider the bi-moment matrix 

hj= I x'y' x''y''dxdy. 

J Jo X + y 

With the change of variables x = rs ,y ^ r{l — s) the integral becomes 

r t^+^+Hl-ty+'+'dt ^ ^{a + ^ + mb + J + l) 
Jo t{^~t) (i + a + j+b + l) ■ ' 

Notice now that the Hankel moment matrix M for the Jacobi polynomials on [0, 1] with weight is 
given by 

/■I 1 

M„- = / x'+^x''+''dx ^ ^ hj =T(a + i + l)M.uT(b + j + 1). (2-2) 

^ Jo i+j +a + b+l ^ ^ ) ^3 \ J I \ J 

This immediately implies the following Proposition. 

m 

Proposition 2.1. Let Pm{x) — Cm.jX^ denote the is the m-th Jacobi orthogonal polynomial then 

i=l 

in n j 

VM :=E^-^r(a + . + l)' ?„(y) ■■-Y.^n, I (2-3) 

are the Cauchy biorthogonal polynomials (not orthonormal) associated with densities x°'e~^ and y^e~y 
and 

Pm{x)qi{y) , dxdy ^ Pm{x)Pi{x)x°'dx = hmSmi (2-4) 



x + y 

where h,n = 27n+a+i ■ result is valid as long as a + 1 > 0, 6+1 > 0, a + b + 1 > 0. 

Remark 2.1. The partition function of the corresponding Cauchy two-matrix model can be easily 
computed using Proposition 2.1 : in [2] it was shown that if c„ are the norms of the monic polynomials 

/I ill • • 11 1 (monic) (monic)\ 

(denoted here provisionally by pn ,qn ) 

p(r°"^H^)9f''"'^(2/)^^^^^da:dy = c„<5„, (2-5) 
2 X + y 

^Here we define the orthogonality on [0, 1] rather than on [—1, 1], The latter is more customary, but the correspondence 
between the two versions is a simple affine transformation of the independent variable x ^ 2x — 1. 
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then Zn = {nl^ n?=o '^r Inspection of the leading coefficients of p„, g„ in Theorem 3.1 shows that 



n\T{a + n + 1) 



T{a + n + l)T{h + n + I) 



r(a + 2n + l) 

n\T{a + n + 1) 
.r(n + 5^)r(n + l + f) 



2n + a + 1 

2 



r(a + n+ l)r(6 + n + l)7r 



(2-6) 
(2-7) 



n-1 



r {^) G{n + l)2G(a + n + ifGi'^fGi'^ + 1)^ G(a + n + 1)G(6 + n + \)n 



r(n+2^) G(a + l)2G(5i±i +„)2G(„+1_^|)2 G(a + 1)G(6 + l)2»("-i)+"° 

(2-8) 

Here (and only in the above formula) G{z) denotes Barnes' G-function, satisfying the relation G(z+ 1) = 
r(z)G(z) , G(l) = 1. 

2.1 The kernels of the correlation functions 

The statistics of eigenvalues of Mi , is expressible in terms of four kernels that can be expressed in 
terms of the CBOPs and auxiliary functions. In keeping with the notation of [3] we introduce the auxiliary 
functions 



Pn{x)x°'e ^Ax 



w-y 



The four kernels to be computed are 



(2-9) 



4ii^,y) 



<\^,y) 



n-1 



^Pjix)q'j^\-y) = 

3=0 

n-1 

^pf\-x)qj{y) = 



3=0 
n-1 



y"'e-y'd y 

+ y + u 

X'°-Q ^ dx (n) , 



K^ix^y) ■■-Y.vfi^~-)qf{-y)-^^ 

3=0 ' ^ 



(„) , .T"^c--'dx' A^^'dy' 1 
Koo'ix ,y) — — — . 

2 x + x' y + y X + y 



(2-10) 
(2-11) 
(2-12) 
(2-13) 
(2-14) 



Each of these kernels can be recovered in general by solving a Riemann-Hilbert problem [3, 2] but in the 
present setting this will be a direct computation. 
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2.2 Facts on Jacobi Polynomials 



Since Jacobi polynomials will be instrumental to our computations, we recall their basic properties. 
Jacobi polynomials are defined by the formulae 



r(a + n + l) 



nl T{a + B + n + 1) ^ \m Via. + ?7i + 1) 



n\ r{a + (3 + n + m + l) f ^-l 



_i Zn + a + p + 1 L (n + a + p + l)n\ 



We shall use the variable 
1 



so that the orthogonality becomes 



z"(l - zfPt''^'> (1 - 2z) Pj°'''^\l - 2z) dz 



1 



T{n + a + l)T{n + (3 + 1) 







2n + a + l3 + l r{n + a + /3 + l)n! 
The case of interest to us is /3 = 0, a = a + b and thus we shall introduce a simplified notation 

r(Q; + n + m + 1) 



P„(z) :=/'("-")(l-2z)=^ 



-^^ m\{n — m)\T{a + m + 1) 



(-2) 



T{a + n + l-u)r(u) du r(n + a + 1) 

Z TT— — — —2J'l 



r{n+l + u)T{a-u + l) 2i-K n!r(a + 1) 
P„,{z)Pn{z)z'' dz ^ 



-n, n + a + 1 
a + l 



(2-15) 

(2-16) 
(2-17) 
(2-18) 



2n-fa + l 

where the integral representation is valid for < |z| < 1. In eq. (2-17) the contour embraces u G IR_ 
(Fig. 1); note that the poles of r(u) for u = —n — 1, —n — 2, . . . are cancelled by the poles of the term 
r(rt + 1 + u) in the denominator and thus the result is a polynomial as a consequence of a simple residue 
computation. 

The Christoffel-Darboux kernel for Jacobi polynomials [10] is de- 
fined as 



Jn{z,w) := ^ 



Pmiz)Pmiw) 



(2-19) 



where hn are the norms squared of P„ in eq. (2-18). The following 
proposition is simple but we could not find it in the literature, and 
since its proof is very close to proofs later on we present it here. 

Proposition 2.2. For /3 = 0, the Christoffel-Darboux kernel for Figure 1: The typical contour for 

Til ■ 1 L J Meiier-type integrals. 

Jacobi polynomials can be expressed as 

n-l 

Kr,{z,w) ^(2j + a + l)Pj{z)P,{w) = 



T{u)T{v)r{n + a + 1 - u)T{n + a + 1 - v) 
2 27ri 27rj I + a — v — u T{n + u)T{n + v)T(a + 1 — u)T{a + 1 — w) ' 



(2-20) 
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valid for < |z|, |w| < 1. Alternatively, introducing the function Gn{z) (a polynomial of degree n — \) 

^' J^2Trir{n + u)T{a + l-u) ^ ' 

we can write ^ 

Kn{z,w)= [ Gn{tz)Gn{twKdt. (2-22) 



At this point we recall the definition of the Meijer G functions, or rather a class of G-functions 
pertinent to this paper. 

Definition 2.1 (Meijer G-functions). Suppose two pairs of natural numbers p < q and < m < q, < 
n < p are given. Then the Meijer G-function is defined by its Mellin-Barnes integral as follows: 



p,q 









-I 




..,bq 







Stt* UUm+i r(i - b, - u) n^^„+i r(a, + u) 



The contour 7, depicted in Fig. 1, is chosen to encircle all the poles of functions T{bj + u) and none 
of the poles of functions T{1 — aj — u) (the implicit assumption is that none of the poles of the former 
coincides with any of the poles of the latter). 

Further properties of G-functions are discussed in the Appendix A. We can identify all the functions 
discussed so far in terms of G-functions. We have 

Proposition 2.3. The Jacobi polynomials P.^ and the function G„ appearing in the representation of 
the Christoffel-Darboux kernel are G-functions with symbols 



Pn{z) = G 



1.1 
2,2 



-a — n,n + 1 



^7i[Z) — ^2,2 



0,-a 

We shall need the following Lemma, whose elementary proof is omitted. 
Lemma 2.1. Let ill = '^''^%Zfr\ttT-''' ■ Then 

Proof. [Proof of Prop. 2.2] From eq. (2-19) and the value of the norms squared given by eq. (2-18) 
we have 

n-1 

^(2j+a + l)P„(z)F„(u;) = 
3=0 

n-1 



du dv z-"w-''r{u)T{v) {2j + a + l)r(a + u)T{a + j + 1 - v) 
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Applying now Lemma 2.1 to eq. (2-24) we obtain 



du dv 



7x7 



27ri 27ri 1 + a — u — v 



du dv 
27ri 27ri 



T{u)T{v)r{n + a + l-u)T{n + a + l-v) 
T{n + u)r{n + v)T{a -u+ l)T{a - v + 1) ^ 

r(u)r(w)r(n + a + l - u)r{n + a + l-v 



dt. 



T{n + u)T{n + v)T{a -u+ l)T{a - i; + 1) 

The term with the 1 vanishes identicaUy because the contours of integration can be retracted to —00 
and z,w € (0, 1), and the denominator does not vanish (a > —1) because the contour around can be 
moved sufficiently to the left so that the denominator has no zero anywhere inside the contour. This 
immediately yields the first formula eq. (2-20). To obtain eq. (2-22) we switch the order of integration 
in the second integral formula above. □ 

It is interesting to consider the scaling limit in which n — > cx) and x, y are rescaled by n^^; using the 
Stirling approximation formula one obtains readily the following Lemma. 



Lemma 2.2. We have the following asymptotic formula 



r(z + 6) 



1 



{5 + p-l){5-p) 



T{z + p) ' y ' 2z 
where the estimate is uniform in each proper subsector. 

Then, using Lemma 2.2 and Proposition 2.2 

-2 „,„-2^ 



0{z- 



|arg(z)| < TT, 



(2-25) 



Kn{xn '^,yn ^) = 
du dv 



x-^y-" 



Tiu)T{v) 



2i7r 2i7T I + a — v — u T(a — u + l)r(Q! — v + 1) 
^ ^ a{a + 1 + xd^ + ydy) ^ q^^--^^ 



a{a + I — u — v) 



2a+2 M a{a + 1 + xd^ + ydy) , ^^^^_2^^ f du dv r{u)r{v)t°'{tx)-^{ty)-''dt 



n 



_,2q+2 



a{a + 1 + xdx + ydy) . _2 

n I .10 



^2 2i'K 2in Jq r{a — u + l)r(Q; — v + I) 
1 



t"Ba{tx)Ba{ty)dt, 



where we have introduced the function Ba defined as 



B^z) := f 



du 

2iTT' 



- z-"/2j„(2yi). 



(2-26) 



(2-27) 



In the above Jq. denotes the Bessel function of the first kind. We also note the integral makes sense for 
any value of z thanks to the super-exponential behavior of the gamma functions. The leading term of 
the kernel in eq. (2-26) is almost the same as the Bessel kernel. Indeed, using eq. (2.2) in [11] we can 
write the Bessel kernel, denoted here by Kb, as 



Ja{Vx)y/yJ'^{^) - Ja{y/y)VxJ'a{Vx) 

2{x ~ y) 



Ja{Vxi)Ja{Vyi)dt. 



(2-28) 



10 



Thus 



A direct inspection shows that 

K^{4x,4y)^] f M2V^t)J^{2^1Ji)dt={xy)^] f B^{tx)B^{ty)dt. (2-29) 

K^{xn-\ yn-^) = n^"^^ ( 1 + + ^ + + ^/^.) + o(,-2)\ ^KsjAx Ay) 

\ n J {xy)2 

This is not surprising since the Bessel kernel arises precisely as the scaling limit of the Jacobi matrix 

model near the hard-edge. 

3 Cauchy-Laguerre Biorthogonal polynomials 

The reason for the use of Laguerre's name is clearly justified by the shape of the weights; nevertheless, 
as we are going to see, the resulting polynomials are more closely related to Jacobi polynomials. Using 
the explicit form of Jacobi polynomials as given by eq. (2-16) and their relation to Cauchy BOPs (see 
eq. (2-3)) we can immediately assert the following 

Theorem 3.1. The polynomials pn of Proposition 2.1 are given by 



T{a + n+ j + l){-zy r(a + n + l) ^ ( -n,a + n + l 

~3)\T{a + j + l)T{a + ] + l " 

du T{a + — u + l)r(u) 



^"^"'^ j!(n-j)!r(a + j + l)r(a + j + l) " n!r(a + l)r(a + 1) ' ^\ a-f l,a-t-l 



(3-1) 
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2i7r r(n + 1 + u)T{a + 1 - u)T{a + l-u) 

where the contour 7 is the same as in Fig. 1 with the poles of T{a + ?i + 1 — w) in its exterior. The 
integral representation is valid for any z ^ 0. The expressions for the polynomials cfn follow from those 
for Pn by exchanging a ^ b. 

Proof. The expression in terms of (generalized) hypergeometric function follows from its definition (see 
16.2 in DLMF); he integral representation is a simple consequence of the Cauchy residue theorem and 

; J_ 
fc! • 

We have two simple corollaries 



the fact that res^r(u) = (-1)'^'^. □ 



Corollary 3.1 (Averages of PmQn)- The average of p^ with respect to the measure x"e ^dx, (g„ with 
respect to y^e^^dy) is 

x''e-^Pn{x)dx = (-1)" = / y"e-^g„(y)dy. (3-2) 







Proof. A simple computation gives: 
Jo f^^ j!(n-j)!r(a + j + l)r(a + j- + l) j!(n - j)!r(a + + 1) ^ 
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It is well known that P,i{l) = (—1)" as a consequence of eq. (2-16). Q.E.D. 
The second corollary is cLii elementary identifica-tion of Pn-i Qn 

respectively, with one of the Meijer 
G-functions. As we will demonstrate the Meijer G-functions play a significant role in the forthcoming 
computations. 

Corollary 3.2. The polynomials Pn and can he identified with the following Meijer G-functions: 



1 



Direct inspection of the leading coefficient of p„, (?„ together with Corollary 3.1 imply the next result 

Theorem 3.2. Let 7r„ := V2n + a + l^ p^g+^'+^^j , ?]„ V2n + a + 1 pg+'^^+^^j - Then the bi- 
orthonormal polynomials Pmin, normalized to have identical positive leading coefficients, are given 
by 

Pn{z) = (-l)"7r„p„, qn{z) = (-l)"77„g„. 
Moreover, 7r„ = a;°e">„(a;)da; and = y^e~'yqn{y)dy. 



3.1 The auxiliary functions pn^ and pn 



Since the definition of qi^^can be obtained from that of pn' by swapping a with b we subsequently will 
focus only on pn'' as defined in eq. (2-9). 

Our goal is to express Pn\z) in terms of the G-functions. To this end we first compute ^ "^^^"^'^^ da; 
for z > 0. 

Lemma 3.1. Let z > 0,a > —l,b > —l,a + b ^ a > —1. Then 



(1) 



Z + X 



dx = e'^G 



2,3 



—n — b, a + n + 1 
0, a, —b 



dx = 



du 



r{a + n-u + l)T{u) 



-dx. 



Proof. By Theorem 3.1 

r°° x^e^^pnix)^^ _ 
Jq z + X 2TTi r(n -f- 1 -f u)T{a + 1 — u)r{a -\- 1 — u) Jq z + x 

where the contour 7 is chosen for this computation in such a way that 5Ru < 1 + a in addition to 5Ru < 1 + a 
required by Theorem 3.1. Furthermore, we move the contour 7 sufficiently to the right to ensure that 
M = a is inside the contour. The change of the order of integration is justified because both single 
variable integrals are uniformly convergent on compact sets of their respective complementary variables 
and the original iterated integral is absolutely convergent. The inner integral is up to a normalization 
the Kummer function of the second kind (see [1] p. 194). A quick way of computing this integral is to 
substitute = e^*^^+^^'*ds then switch the order of iterated integrals to get: 

/•oo a—Up—x 

dx = e^r(a + 1 - u)z''-"r(u - a; z), 

Z + X 
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where T{u — a; z) is the incomplete gamma function. To conclude this part of the computation one only 
needs to use T{u — a; z) — T(u — a) + "f{u — a; z) — T(u — a) — 1-^1(1* — a,u — a + 1; — z) which 

expresses the complementary incomplete gamma function 7 in terms of the hypergeometric function iFi 
(see [5] p. 266, formula (22)) thus obtaining 

/ — dx = e'^r(a + 1 - u)T{u - a)z""" - c'^r(a - u) iFi{u - a,u - a + 1; -z). 

Jq z -\- X 

Putting the first term on the right side back into the contour integral gives the claim if one uses the shift 
formula eq. (A. 2). To finish the proof we need to show that the contour integral of the second term 
vanishes. To this end we write explicitly the resulting contour integral obtaining, after omitting e^ and 
performing elementary computations with the gamma functions, 

du T{a + n — u + l)T{u) 1 / -.j 

^ 2TTi r{n + 1 + u)r{a + 1 - u) ^ (u - a + ^ ^' ' 

The series converges uniformly on 7 hence to prove that this integral is zero it suffices to prove 

du T{a + n - u + l)r{u) 1 _^ -Ql 

2m r(n + 1 + u)T{a + I - u) (u - a + j) ^ ' ^ ' ' ' ' ' 

The integrand is a meromorphic function with poles at u = {0, —1, • • • , ~n\ U {a — j} which, accord- 
ing to our choice of the contour 7, are all inside the contour. We can now extend this contour by 
adding a large circle with the center at u = and radius r and observe that the integrand on such 
a large circle is, in view Lemma 2.2, 0{u~^). This implies that in the limit oi r ^ 00 the integral 
Xy |^ r[i+t+»)rta+i-l) (u^i+j) = [ sum of residues of the integrand outside of 7] = 0. 

□ 

The formula for p"n\z) follows now easily from Lemma 3.1 and Theorem 3.2. 

Theorem 3.3. For z G C \ M+ the auxiliary functions of the first kind pn ^ have the following Meijer 
G-function representation 

^(,)^^^ J- pr.ix)xy-^dx ^ ^-.(_,)„+i^^g2,i -n -b,a + n + 1 _ \^ 

A similar expression holds for by interchanging a ^ b and iTn <^ r]n. 

To see the relation between p^ ^ and Pn wc formulate the following equivalent representation of pn ^ 
whose gamma part is identical to that for p„; see eq. (3-1) 

Corollary 3.3. For z G C \ M+ the auxiliary functions of the first kind admit an equivalent integral 
representation: 

- r_1^n+l^ f'^ r(l + n + Q;-M)r(M) 



Pniz)-e'i 1)" l2mr{l + a-u)r{n+l + u)r{l + a-u)^ \in(7r(u - a)) ^^"^^ 
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14 / -n+a,n+l \ 
2,3 \^ Q-a,~a ' 



Proof. The G-function occurring in the representation of |5„ is G2'3 ^ J*^"'"^^ —G 
Likewise, in view of Theorem 3.3 and the shift formula eq. (A. 2) the G-function in the representation 
of p"n^ can be written as {—z)°-g1^\ ~"~o'"a^ ^ Switching the first and the second coefficient 
in the bottom line of the latter G symbol can be done if one uses Euler's reflection formula to write 
T(u — a) = T^TFT-, , / which implies the claim. □ 

^ ' r(l+a.— ti) sin(7r(u— a)) ^ 



Remark 3.1. It is instructive to verify that pn\z)-^- ~pn'{z)^ = — 27rie ^z°'Pn{z). Indeed, let z > 
then 



(-z)r"-(-^)° 



gZ7r(ii— a) ^—i7r{u—a) 



z''-"2isin(7r(u-a)), 



since the cut for z^ is by convention (see A) along the negative real axis of z. Thus by the corollary above 
we get: 



pl,^>iz)+~p^^^>{z) 



du 



Tin + a + l-u)T{u) 



^2T:i T{a + 1 - u)T{n + 1 + u)T{a + I - u) 
which by Theorem 3.2 implies pn\z)^ -~pn\z)- — ~2'Kiz°'e^^Pn{z). 



z'^"" 2i siri(ji4*r==-a)X- 



^iri43ii(«-==Tff)' ' 



(3-4) 



Theorem 3.4. For z € C \ M_ the auxiliary functions of the second kind admit the G-funetion represen- 
tation 



oo poo 



— n, a + n + 1 
0, a, b 



(3-5) 



Proof. Let us express pl^"^ in terms of pil^ as 



— \ — p«'(-y)dy- 

z + y 



Since (-y) = (-1)"+V„ e^ G^ ( -"^l'^^-^' | v) 
p(2)(z)^(-l)"+V„ 



we obtain 



V ^2,1 ( -n-h,a + n+l 



/o z + y 
by eq. (A. 4) and eq. (A. 2). 



a, 0, —6 



2/ dy=(-l)"+V„ G 



2,3 



-n, a + 71 + 1 
0, a, 6 



□ 



Again to see the relation between pn^ and pn^ we establish an equivalent representation 



Corollary 3.4. For z G C \ M_ f/ie auxiliary functions of the second kind admit an equivalent integral 
representation: 



J2)/^ / iN„+i f r(u)r(M + a)r(l + n + &-M) _^fc_„ 



2ni r(l 6 - u)r(a + n + 1 + u) sin(7r(u - 6)) 
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— n — b,a+n+l 
-b,a,0 



Proof. Using the shift formula given by eq. (A.2) we write GI'I ( ""o"^^^^ z) = z''GI[1 ( 

On the other hand, in view of Theorem 3.3, the G function in the representation ofpi^^ is Gj'l^ ^ ""o^a^-b"^^ 
Switching the first and the last coefficient in the bottom line of the G symbol can be done if one uses 
Euler's reflection formula to write T(u — b) = wttt — n^- , — tt which implies the claim. □ 



Remark 3.2. By analogy with the case we easily verify the jump of pn' (z) for z < to be 

p^n\z)+ —pn\z)^ = 2'iTi{~z)^e'^pn\z). Indeed, let z < then 



J.'K{h — U) 



-iixih—u) 



= (-z)''-"2isin(7r(6-M)), 



since the cut for z'^ is by convention (see A) along the negative real axis of z. Thus, again, by the corollary 
above we get: 



I f —n — a + 71 + 1 



p^^\z)^^p^^\z)^ = (-z)''(-l)"+V„G2,3 y 
which by Theorem 3.3 implies pn\z)j^ —pn\z)^ = 2ni{—z)''e''pn\z). 

3.2 The kernels for finite n 

In [3] we proved a variety of generalizations of Christoffel-Darboux identities. In this section, however, 
we show that in the case at hand it is possible to compute the associated kernels directly and in an 
elementary way from the formulae above. We will begin with i^go'' ^ defined in eq. (2-10). 



Theorem 3.5. The principal kernel KQQ{x,y) is given by 



{tx)Gb.,n{ty), 



where 





du 



.T-" r{u)r{n + a + 1 - u) 
2iTrT{c+l-u) r{n + u)T{a - it + 1) 



= G 



1,1 

2,3 



—a — n,n 
0, — c, —a 



(3-6) 
(3-7) 



Proof. We start with writing the formula for the kernel in terms of G-functions. By Theorem 3.2 

while Christoffel-Darboux kernel 



K, 



(n) 
00 



i^.y) - E;=oi2j+»+^)Gl:l{-o'::i'-i' I-) G^^{-^:^^ \y) 



for Jacobi polynomials is given by Kn{x, y) — X^^Co (2j + ct + 1)G2'2 "o -q^^ 
(see Propositions 2.2 and 2.3). This means that the integral representations obtained for Kn(x,y) are 
only modified by extra factors Y'{i+a~u) ' r{i+b-v) coming from the additional third index appearing in 
the bottom of the respective G-symbols oi pj{x),qn{y) which proves the claim. □ 



G 



1,1 / -a-].j + l 

0,-Q 
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Remark 3.3. We can write explicitly the double integral representation of the kernel by copying the 
analogous formula from Theorem 2.2 and attaching additional factors Y(i+a-u) ' r{i+b-v) explained in 
the proof above. The formula reads: 



du dv 



T{u)T{v)T{n + Q + 1 - u)r{n + a + 1 - v) 



(3-8) 



2m 27ri T{a+l- u)T{b + 1 - v){l + a - v - u) T{n + u)T{n + v)T{a -u + l)T{a - v + I)' 

Note that this integral representation, as opposed to the one appearing in 2.2, is valid for any x,y ^ 0. 
Observe also that the function Gc,n{x) — Gj'^^ ^ o^"c -a ^) ^'^ natural generalization of the function 

Gn{x) — G2 



1,1 / — Q — r 



2,2 I 0,-Q 



xj known from Theorem 2.2. 
We are now ready to compute the remaining kernels. Using Lemma 2.1 we find 
Theorem 3.6. The kernel K^^^x^y) is given by 

K^y;\x,y)=x''e'= eGaAtx)Gb,n{ty), where 



Gc,n{x) 



du r(u - c)r{u)r {n + a + l - u) _„ _ „2.i 

x — Gr 



27rj 



r(a + l-M)r(n + u) 



2,3 



—a — n, n 
0, — c, —a 



Likewise, by symmetry, the kernel K'^{x,y) is given by: 



K^^\x,y) = y'ey I rG,.„(te)Gfc,„(<2/) 
Jo 



with Gc,n defined in eg. (3-10) and Gc.n defined in eq. (3-7). 

Proof. From the definition of k[q^ in eq. (2-12), Theorem 3.3 and eq. (A. 2) we obtain 

n-l 



(3-9) 
(3-10) 



(3-11) 
(3-12) 



K[^\x,y) ■.^-Y^pf\-x)q,iy) = 



j=o 



n-l 



x°e-5](2j + a + l)G^:,^ 
3=0 



0, —a, —a 



G 



1,1 
2,3 



0,-a,-b 



We observe that the dependence on j is identical to the one covered in Lemma 2.1. Hence by that lemma 

r-l 



K['(,\x,y)=x'^e- 



ft-Gli 
Jo 



t"G: 



2,1 
2,3 







—a — n,n 
0, —a, —a 



tx G 



1,1 

2,3 



—a — n, n 
0, —a, —b 



ty\dt- 



tx ] G°i 



-6 



tx Gi 



—a — n, n 
0, —a, —b 



since the second term, or more precisely Gq'^ ( _^ ty\, vanishes by eq. (A. 5) 



ty ]dt, 

□ 
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Likewise, 

Theorem 3.7. The Kernel K[^\x,y) is given by 

,(n) 



x + y 



(3-13) 



with Gc.n defined in eg. (3-10). 



Proof. We plug the expressions for pn \ qrP of Theorem 3.3 into the definition of KYi' (2-14) to get 

n— 1 



(") 



cVe^+^^(2j + a + l)G^:J 

3=0 



-a- i,j + 1 
0, —a, —a 



G: 



2,1 
2,3 



0,-5, -a 



Again, the j-dependence of this expression is covered by Lemma 2.1 and applying it we obtain: 



1 



KYl\x,y) + ^^x'^y'e-+y{ i"G 



x + y 



21 I —a. — 71, n 



'2,3 



0, —a, —a 



tx Gi 



2,1 / -a-n,n 
0, —6, —a 



tyUt- 



tx 1 Gn "i* 



-6 



%)dt}. 



By inspection GJ'J ( _^ tx) = (te) °e hence an elementary computation gives the final answer 



1 (e-(^+2/) _ 1) 

2;ayh X + y 

which implies the claim. 



-a;°/e"+^{ / i"G: 



'2,3 



a^2,i ( -a-n,n 
0, —a, —a 



ia; Gi 



2,1 / -a - n,n 
0, — &, —a 



□ 



Remark 3.4. Similar to the integral formula for the kernel Kqq (.see eq. (3-8)) all the remaining 
kernels have analogous integral representations in terms of the double path integral, all involving the 
kernel 



I + OC — U — V ' 

3.3 The Meijer-G random point field 

In order to describe the statistics of the Cauchy matrix model near a: = = y we recall the definition of 
the kernels for the correlation functions, that is, the K and H kernels defined by eq. (1-6), eq. (1-5)) 
respectively. 
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3.3.1 Asymptotics of the kernels K^Z^ 

From the explicit integral expressions for the functions Gc,n,Gc^n and the kernels in Theorems 3.5, 3.6, 
3.7, we shall now derive the behaviour under the rescaling 



X 



c 



y 



i2' 



n — oo. 



Remark 3.5. It should he mentioned here that if we had started with a model in the form 

du Ml , M2 ) = — dMi dMa , . . , ^ 

' ^ Z„ dct(Mi+M2)" 



(3-14) 



(3-15) 



then the relevant rescaling would have been x i— > xn ^ , because the n scaling in the exponential is simply 
absorbed by a rescaling x — nx. 

Theorem 3.8 (Kernels in the asymptotic regime). The following asymptotic estimates for n — > cx) hold 
uniformly over compact subsets of the independent variables 



= :Gc(C) 



n-"-^Ge,„(C«-') = 



du 



T{u) 



2m r(c + 1 - u)T{a + l- u 

l^ ;:i";f^^ r"+c^(n-^), 

ZTTl 1 (a + 1 — u) 



■r"+o(n-i). 



(3-16) 
(3-17) 



= :Ge(C) 



or in terms of G-functions 



n—'G\'l 



—a — n, n 
0, — c, —a 
~a — n, n 
0, — c, —a 



0,3 



0, — c, —a 



G, 



0,3 



0, — c, —a 

Correspondingly, the asymptotic behavior of the four kernels is 



Cj -)-O(n-i), 



'-K'->{Cn-\in-^) 



Ga{tC)Gbmt"dt, 



n 



-2a-2 7^(") 



K'^i\Cn'\in-') ^ e / GM)GM)t^At, 
Jo 

i-2''-2if(^)(Cn-2,C«-2)^C / Ga{tC)Gtmt"dt, 



n-^KrACn-^^n-^)^C^'' / G,(tC)Gf,(te)i"dt 
and the convergence of the above limits is all within 0{n~^). 



(3-18) 
(3-19) 
(3-20) 
(3-21) 
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Proof. The proof of this theorem is immediate from the Stirhng approximation formula in Lemma 2.2, 
the exphcit integral expressions for Gc,mGc,n in eq. (3-7), eq. (3-10) and the convolution form for the 
kernels in Theorems 3.5, 3.6, 3.7. □ 

Now we turn our attention to the H kernels defined in eq. (1-5). Using the rescaling (3-14) and the 
asymptotics of the kernels in Theorem 3.8 we find the following expressions 

(3-22) 
(3-23) 
(3-24) 
(3-25) 



lim n 

n— >C30 


-2 Tjin, 






lim n" 


-1 rr("; 


\x,y) 


= Q^e 


lim n" 

n— >oo 








lim n~ 

n— ^C30 









Remark 3.6. The scaling factor n in front of the kernels Hj^l^ is natural since kernels transform as 
A/dxdy. 

The limiting kernels above define a novel class of random point fields. 

Proposition 3.1 (Meijer-G random point field and universal class). In the scaling limit the correlations 
of the eigenvalues of Mi , M2 are determined by the two-level random point field with kernels 



eoo(C,0 = CT/ Ga{tOGbmt"dt, ^oi(C,0=CT/ GaitOGbmf'dt, (3-26a) 
Jo Jo 

GioiCO = r Ga{tOGM)t"dt, GiiiCO = f' GaitOGbmf'dt^-^ (3-26b) 
Jo Jo S + s 

where S,,C ^ 

Remark 3.7. It is a nontrivial fact that the correlation functions defined using the kernels in Proposition 
3. 1 are all positive and define an actual random point field: of course this is hard to see directly from the 
definition of the kernels; it nevertheless follows from the fact that this is the limit of a well-defined point 
field. 

If we absorb the powers in the definition of Ha, Ha as in Definition 1.1 then the formulas becomes 
even more symmetric. 

Corollary 3.5. In the scaling limit the correlations of the eigenvalues of Mi,M2 are determined by the 
two-level random point field with kernels 

Goo{C,0= [ HaitOmmdt, goi{C,0= f Ha{tC)Ht{tOdt, (3-27a) 

Jo Jo 

GioiCO^ C Ha{tOHM)dt, Gii{C,0 = f' Ha{tOHt{tOdt - (3-27b) 

Jo Jo + s 
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where ^, C G ^+ ^.i^-d the kernels are given by the Meijer G functions 

' \c, 0, — a + c / \c, 0,— a + c 

Thus we have completed the proof of Theorem 1.2 stated in the introduction. 

Remark 3.8. It is easy to see by direct computation that the functions He, He are connected to hyperge- 
ometric functions 0^2 ■ For example 



2° 0^2 (i+a,l+a'-^) 



- r(i + «)r(i + a) ' ^'-''^ 

while Ha{z) is, for a ^ 0, 1, • • • , a linear combination of two hypergeometric functions of type 0^2 

r (a) 0F2 ( , ^) r (-a) z- 0F2 ( , z) 

^ r(i + 6) r(i + a) ^ ' 

The remaining cases a = n — 0,1, ■ ■ ■ are handled by taking limits a ^ n. For example 
3.3.2 Expressions of the kernels in "integrable" form 

We note that the functions Ha{t), Ha{t) solve differential equations of the third order; this fact is implied 
by their identification as Meijer-G functions. A simple derivation is as follows: denote by A{u) the 
rational expression of F functions in the integrand of He or H^, then a direct inspection shows 

Aiu + l) = Tu{u + c)ia-c-u)A{u) ^ Ac(Ac + « - c)(Ac - c)/(C) = tC/(C), (3-30) 

where = C^i the upper sign is for He, the lower for He respectively. Now, let c G {a, b} and denote 
by /t(C) f{tC), 9t{0 -=5(^0: s-iid / be one of Ha or Ha and g one of Hb or Hb. Since A is scale 
invariant we have 

5t(0Ac(Ac + 6)(A^ - a)/4(C) - TtCMOgtiO, (3-31a) 
/t(C)Ae(A5 + a)(A^ - b)gtiO - Ti^/t(C).9t(^)- (3-31b) 

Adding these two equations together, dividing by t and using that Ac^ft{C) — At/t(C) (and similarly for 
9t{0) we find 

dt [gtA^Jt - AtgtAJt + f^A^.g^ + (6 - a)[gtAtf - ftA^gt) - abftg^] = tC + T0/t5f (3-32) 
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Integrating with respect to t from to < r we obtain 



TC + tO [ ftgtdt = ftA^.gt ~ AtftAtgt + gtA^ft + {h-a) {gtAtft - ftAtg) - ah ftgt 

[f(tO,Acf{tO,Alf{tC)] 



0+ 



—ab a — b 1 
b-a -1 
1 



Ae.9(te) 



A?.9(te) 



, (3-33) 



t=o+ 



where in the last hne we used again the scale invariance of A. The last computation motivates the 
following definition 

Definition 3.1. Let f — /(C) and g = g{£,)- Then the point-split bilinear concomitant is defined as 
B{f, g) := fAlg - {A^f){A^g) + gA^ +{b-a) {gA^f - fA.g) ~ abfg = 



—ab a — b 1 
b-a -1 
1 



.9 

A^g 
A|5 



(3-34) 



If C = suitably chosen to make the left side of eq. (3-33) vanish, then eq. (3-33) implies 
i9{i3(/, g) — 0. This is a form of the bilinear concomitant [7] for equations in duality, which explains our 
naming convention. 

From this point onward we will be interested in r = 1. The evaluation of the right hand side of eq. 
(3-33) a,t t = 0+ is case dependent. 



Case Ha, Hb. If / = Ha, g = Hb then the evaluation at i = vanishes: indeed by eq. (3-28) 

1 



r(i + c)r(i + a) 



o(C) 



(3-35) 



and it suffices therefore to verify S(C°,C ) — and 



lim m(tC)''+ , (i£)^ ) = lim e+^^+'BiC^^Cn = 0, for 1< fc + / , 

since a + 1 > 0. 

Case Ha, Hi, or vice versa. This case is similar to the previous case; one uses eq. (3-29) and analyzes 
cases. The relevant points to remember are that by our assumptions a + l>0, 6+l>0 and a + 1 > 0. 
In all cases B vanishes at t = 0+. 
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Case Ha, Hi,. In this case the \imt-fO+ B{Ha, Hb){t(,t^) = —1. This follows from the fact that in the 
expansion of Ha (or Hb) there is one more pairing that the previous cases did not have. This is the 
pairing of the leading term in eq. (3-29), namely in '"("^ °^T{i+bj^ '^ ^'^ ' with its analog in H^. Extracting 
the leading term leads to i3(i^, rg^) = -a ^(i7a)ra+fc) = 

It is now elementary to substitute the above results into eq. (3-33) and, further, into the definition 
of G kernels given by eq. (3-27a). 



Proposition 3.2. Let B be the point-split bilinear concomitant defined by 3.1. Then the correlation 
kernels of the Meijer-G random point field satisfy 



B{Ha{0,Hb{0) 

B{Ha{C),Hb{0) 



Gii{CO = 



g(g.(C),gb(0) 

B{Ha{C),Hbm 



(3-36a) 
(3-36b) 



A The Meijer-G functions 

For the convenience of the reader we recall the definition of the Meijer-G function 



p,q 



f ai,. 


. . , Op 




-I 


\ bi,. 


■■,bq 







2'^^ n =™+ir(l-&. -")nU+ir(a 



(A.l) 



The contour 7 is a contour that leaves the poles of the r{bj + u)'s to the left and the poles of the 



r(i 



u)'s to the right (the implicit assumption is that none of the poles of the former coincides 



with any of the poles of the latter). The contour extends to infinity in with | arg(w)| < ^ (in the right 
half plane) ii p > q, while it extends to infinity with | < arg(u) < ^ (in the left half plane) if p < g. It 
is the latter case that is relevant for us. We assume therefore that the branch of z~" is chosen in such a 
way that the cut is for z < 0. There is an extensive literature on G- functions (see e.g. [8, 9]). We collect 
several useful formulas used in the paper. 



1. Shifting formula: 



p.q 



^ m.n 

p,q 



bq -I- CT 



where +a means shifting every entry by a. 
2. Integrals containing products of two G- functions: 



(A.2) 



Since G^q 



-bi,. 
-ai, . 



is proportional to z'^ i-fo(l + b — a, — z) = 



(z + y)' 



_ m.n 



iJ-y 



Via) 



G 



771+1, n+l 



1 — a, a. 



-bm+i, ■ ■ ■ , —bi 

~^n+l 1 ■ ■ ■ ; ^] 

we obtain [8] 



rjU! 
(A.3) 



(A.4) 
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3. For p < q 



0. 



(A.5) 



B Correlation functions 



In this appendix we recall the definitions and formulas used to study the correlation functions for generic 
Cauchy matrix models (no restrictions on the measures). We use the notation that is slightly different 
than the one used in [2] to accommodate the needs of the present paper. Given two measures da{x) = 
a{x)dx, d/3{y) = /3{y) with densities a{x),/3{y) respectively, both supported on the positive half-line 
the Cauchy matrix model is the probability measure on the pairs (Mi, M2) of positive semi-definite NxN 
matrices given by 

a(Mi)/3(M2)dMidM2 



dAi(Mi,Af2) 



Zn det(Mi Ma)^ 



where a{M) (or 13(B)) stands for the induced measure on the spectrum of M, i.e. nja(a::j) where Xj 
are the (positive) eigenvalues of M. Let us introduce two families of polynomials {pj{x),j = 0, 1, ..} and 
{qj{y),j = 0, 1, ..} which are biorthonormal with respect to the pairing with the Cauchy kernel 

p^{x)a(x)qi,(ij)^{y) 



B2 



-AxAy = bj^k, 



x + y 

with the technical proviso that the leading coefiicients are identical and positive to render all polynomials 
unique. Following [2], but changing slightly the notation, we introduce four kernels used in the present 
paper 



n-l 



K^;^(x,y) :^ 1 4^\x,y')^-^, Ki^\x,y) :^ [f 4o\^\y') f2%(K^''^ " ^ ' 
J y + y J J {x + x'){ii + y') x + y 

With this notation in place the correlation functions for r eigenvalues xi^. . . ^Xr of Mi and s eigen- 
values 2/1 , . . . , 2/s of M2 can be shown to be given by 

r s 

TZ^'''''\xi, . . . ,Xr;yi, . ■ . ,ys) = ]^a(a:j)dxj J]^ /3(yfe)d2/fc det 

j=i k=i 

It is easy to write the correlation function as one determinant by including the products in front of the 
determinant above as determinants of diagonal matrices. This is not a unique procedure but the choice 
that works naturally for the present paper is the following. We write symbolically 

a(x) Ol rKoi(x,x) Koo(x,y)l \l 
I Kii(y,x) Kio(y,y) /3(y) 





[KQg{xi, J/j)] ,,<s 




[K[o\y^,yj)],^^<, . 



TZ'^'''''\xi, ...,Xr;yi,.-.,ys) = det 
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where a{x) = diag(a(a;i), • • • , a{xr)), diag(^(yi), • • • , I3{ys)) respectively. Hence 

n^'''^\xi,...,Xr;yi,..-,ys) = det 

This leads to the definition of new kernels: 

Hi,^\x,y) := a{x)l3{y)K^^\x,y), H^i\x,y) := a{x)K^f {x,y), (B.2a) 

Hi^\x,y) := l3{y)Ki^\x,y), H^^\x,y) := K^\x,y). (B.2b) 



a(x)Koi(x,x) 
Kii(y,x) 



a(x)Koo(x,y)/?(y) 

Kio(y,y)/3(y) 
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